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We derive a general expression for the probability of observing deviations from statistical isotropy 
in the cosmic microwave background (CMB) if the primordial fluctuations are non-Gaussian and 
extend to superhorizon scales. The primary motivation is to properly characterize the monopole and 
dipole modulations of the primordial power spectrum that are generated by the coupling between 
superhorizon and subhorizon perturbations. Unlike previous proposals for generating the hemi¬ 
spherical power asymmetry, we do not assume that the power asymmetry results from a single large 
superhorizon mode. Instead, we extrapolate the observed power spectrum to superhorizon scales 
and compute the power asymmetry that would result from a specific realization of non-Gaussian 
perturbations on scales larger than the observable universe. Our study encompasses many of the 
scenarios that have been put forward as possible explanations for the GMB hemispherical power 
asymmetry. We confirm our analytic predictions for the probability of a given power asymmetry 
by comparing them to numerical realizations of CMB maps. We find that non-local models of non- 
Gaussianity and scale-dependent local non-Gaussianity produce scale-dependent modulations of the 
power spectrum, thereby potentially producing both a monopolar and a dipolar power modulation 
on large scales. We then provide simple examples of finding the posterior distributions for the 
parameters of the bispectrum from the observed monopole and dipole modulations. 


I. INTRODUCTION 

It is very tempting to try to use large-scale features 
of the primordial fluctuations, such as the hemispherical 
power asymmetry, as a clue toward primordial physics. 
Such signals are both intriguing (maybe they say some¬ 
thing about the beginning of inflation) and statistically 
unfortunate (i.e., not all that unlikely to be a feature 
of a particular realization of Gaussian, isotropic fluctua¬ 
tions). Many studies of the observed power asymmetry 
and its statistical significance have been reported using 
the WMAP data [IH3] and the Planck data [H-^ . A va¬ 
riety of possible explanations for this asymmetry have 
been discussed (see e.g. [lO]): several of the most in¬ 
triguing ideas use superhorizon fluctuations to generate 
the asymmetry, either by using non-Gaussianity to couple 
them to observable perturbations [nH3i] or by postulat¬ 
ing some different primordial physics that precedes the 
usual slow-roll inflation |551 - I55| . Alternatively, one can 
postulate scenarios that are fundamentally anisotropic 
on the largest scales 

It is quite general that if the primordial fluctuations 
are non-Gaussian, the likelihood of observing statistical 
anisotropies in our cosmic microwave background (CMB) 
changes. Although one might expect that isotropy and 
Gaussianity are independent criteria for the statistics of 
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the primordial fluctuations, this distinction is not actu¬ 
ally clear when we only have access to a finite volume of 
the universe [UHIS]- If tde primordial fluctuations are 
non-Gaussian, the observed large-scale “discrepancies” 
from the simplest isotropic, power law, Gaussian fluc¬ 
tuations need not be a signal of a special scale, time, or 
feature during the primordial (inflationary) era. Instead, 
they may simply be a consequence of cosmic variance in 
a universe larger than the volume we currently observe, 
filled with non-Gaussian fluctuations. 

In this paper, we present a single framework to calcu¬ 
late the distribution of expected deviations from isotropy 
in our observed sky from any model with non-Gaussian 
primordial fluctuations. This framework incorporates 
most successful proposals for generating the power asym¬ 
metry, even some that were not originally formulated as 
non-Gaussian models. The reason is that if the assump¬ 
tion of statistical isotropy is maintained, any explanation 
of the temperature power asymmetry can be modeled by 
assuming a fluctuation in a long wavelength modulat¬ 
ing field that couples to some cosmological parameter 
relevant for determining the CMB power spectrum (the 
fundamental constants, the scalar spectral index and the 
inflaton decay rate, for example uni HU mi 133 SSI 117]). 
Since the asymmetry is well modeled by a (small) range 
of shorter scales all coupling to a longer scale this is, by 
definition, a sort of non-Gaussianity. 

Using our analytic framework together with numeri¬ 
cal realizations, we show that isotropy violations in our 
observable universe are more likely in models where sub¬ 
horizon modes are coupled to longer wavelength modes. 
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In some models they are so much more likely that the 
extent to which our observed sky is isotropic can be in¬ 
terpreted as a constraint on non-Gaussianity. 

In addition, we find that scale dependence of such mod¬ 
ulations is a generic outcome of non-Gaussianity beyond 
the standard local ansatz. This result is important be¬ 
cause a major hurdle in model building for the power 
anomalies has been their scale dependence, as the ob¬ 
servations suggest a very sharp decrease in the power 
asymmetry amplitude at smaller scales [3 |48l |49] . Scale 
dependence of the type observed occurs whenever modes 
closer to the present Hubble scale couple more strongly 
to superhorizon modes than very short wavelength modes 
do. This can be the case even for scale-invariant bispec¬ 
tra. For example, since equilateral-type non-Gaussianity 
peaks when the momenta configurations are nearly equal, 
the modulation from an equilateral-type bispectrum is 
biggest at large scales and quickly dies off at smaller CMB 
scales. 

Further, we will see that any non-Gaussianity used 
to explain the scale-dependent power asymmetry will 
also produce a scale-dependent modulation of the power- 
spectrum monopole. This means that the observed low 
power on large scales and the dipole power asymmetry 
can both provide evidence in favor of non-Gaussian fluc¬ 
tuations. To quantify the significance of these signals, 
we perform a parameter estimation of the non-Gaussian 
amplitude and shape using both monopole and dipole 
modulations, and estimate Bayesian evidences against 
the Gaussian model. 

Finally, we will show that scale-dependent non- 
Gaussianity eliminates the need for enhanced inhomo¬ 
geneity on superhorizon scales to generate the observed 
asymmetry. The most common method for introducing a 
dipolar power modulation is to postulate the existence of 
a large-amplitude superhorizon fluctuation in a specta¬ 
tor field during inflation that then alters the power spec¬ 
trum on smaller scales via local-type non-Gaussianity. 
(A multiple-field model of inflation, or one that other¬ 
wise breaks the usual consistency relation, is required 
because superhorizon perturbations in the standard in- 
flaton field cannot generate an asymmetry [^.) Dubbed 
the Erickcek-Kamionkowski-Garroll (EKC) mechanism 
this approach has been expanded upon and refined 
several times since its inception m [HHini [1H31]. in 
these analyses, WMAP and Planck bounds on local-type 
non-Gaussianity forced the amplitude of the superhori¬ 
zon perturbation to be much larger than predicted by 
an extrapolation of the observed primordial power spec¬ 
trum to larger scales. Although possible origins for this 
large-amplitude fluctuation have been proposed, such as 
a supercurvature perturbation in an open universe [19] , a 
bounce prior to inflation |31j . a non-vacuum initial state 
pS] , topological defects [Mils, or deviations from slow- 
roll inflation |20l[22]) it is largely taken as an ad hoc addi¬ 
tion to the inflationary landscape. The scale dependence 
of the asymmetry often requires an additional elabora¬ 
tion to the theory, either in the form of scale-dependent 


non-Gaussianity [HI |T31 [TSl US or isocurvature 

fluctuations [HEolEl]. 

Dropping the idea of a special large superhorizon fluc¬ 
tuation and instead starting with scale-dependent non- 
Gaussianity changes this picture. Importantly, the con¬ 
straints on the amplitude of the non-Gaussianity on large 
scales is rather weak; the Planck bound of = 

2.5 ± 5.7 at 68% CL [5T] assumes a scale-invariant bis¬ 
pectrum. To generate the observed power asymmetry, 
we only need non-Gaussianity on the scales that are 
asymmetric. For example, if we consider CMB mul¬ 
tipoles up to £ ~ 100, the WMAP 5 year data give 
/nl = —100 ± 100 [HI- (This constraint has not appre¬ 
ciably changed since WMAP 5: see also Figure 11 of m 
for the most recent plot of Planck’s fmax-dependent con¬ 
straints.) Many arbitrary choices for £niax in the range 
[40, 600] (and different binning schemes to study the scale 
dependence) can be found in the literature for comput¬ 
ing the large-scale power asymmetry amplitude and its 
significance [3 [1130. We will use = 100 as our 
fiducial value to compute the asymmetry for our numer¬ 
ical tests and statistical analysis later. This also allows 
us to use the WMAP5 large-scale bispectrum constraints 
computed with fmax = 100. 

The weaker constraint on /nl implies that it is possible 
to generate the observed asymmetry without enhancing 
the amplitude of superhorizon fluctuations, as was noted 
by [13, which considered the power-spectrum modula¬ 
tions generated by anisotropic bispectra. We note, how¬ 
ever, that the dipolar bispectrum considered there does 
not generate a power asymmetry because the k ^ —k 
symmetry of the power spectrum forbids any modulation 
to the power spectrum from bispectra that depend on odd 
powers of the angle between the long and short modes; 
this point was clarified in [53j . In this work, we will focus 
on isotropic, scale-dependent non-Gaussianity (although 
our framework can easily be extended to include funda¬ 
mentally anisotropic models), and we will show that this 
is sufficient to generate the observed asymmetry without 
enhancing the amplitude of superhorizon fluctuations. 

The plan of the paper is as follows. In the next section, 
we use the usual local ansatz to demonstrate the validity 
of our analytic calculations of the statistical anisotropies 
expected in non-Gaussian scenarios and to illustrate sev¬ 
eral of the key conceptual points relating non-Gaussianity 
and anisotropies. We also show how our framework en¬ 
compasses the EKC mechanism and demonstrate that ex¬ 
otic superhorizon perturbations are not required to gen¬ 
erate the observed power asymmetry. In Section m we 
test our analytic calculations for monopole and dipole 
power modulations using numerical realizations of CMB 
maps. We also introduce and discuss our parameter es¬ 
timation and model comparison methods. In Section jTV] 
we investigate non-Gaussianity beyond the local ansatz 
and in particular, consider a representative model that 
generates features that are fully consistent with con¬ 
straints on the isotropic power spectrum and bispectrum. 
We discuss and summarize important aspects of our work 
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and conclude in Section]^ The appendixes contain tech¬ 
nical details. 


II. ILLUSTRATING THE CONNECTION 
BETWEEN NON-GAUSSIANITY AND 
ISOTROPY 

We assume that at some early time (after reheating 
but prior to the release of the cosmic microwave back¬ 
ground radiation) a large volume of the Universe (Vl) 
contains adiabatic fluctuations described by isotropic but 
non-Gaussian statistics. To compare predictions of this 
model with observations, we are interested in the statis¬ 
tics of the fluctuations in smaller volumes, Vs Vl, 
that correspond in size to our presently observable Hub¬ 
ble volume. We will first consider the usual local model 
with constant /nl for simplicity; we will present more 
general results in a later section. 


A. The local model 


Suppose the Bardeen potential <i> is a non-Gaussian 
field described by the local model: 

$(x) = + /nl - (</'(x)^)) , (1) 

where /)(x) is a Gaussian random field. When the large 
volume is only weakly non-Gaussian, the power spectrum 
observed in our sky, P$ 5 (fc,x), will be related to the 
mean power spectrum in the large volume, Pc/,{k), by 


P<l,s(A:,x) = P^{k) 


1 + 4/nl 


(27r)3 


4>(ke)e 




( 2 ) 


where the radial integration for ki is confined to 
|fcf| < 7r/rcmb when the GMB spectrum is the quantity 
of interest. Our conventions for the power spectrum are 
stated in Appendix and Appendix provides the 
derivation of this equation. Any particular model for 
generating the fluctuations in volume Vl should provide 
a well-motivated lower bound on the ki integral (e.g., 
from the duration of inflation). However, not all shifts 
to local statistics are sensitive to the full range of the 
integral; for local-type non-Gaussianity, as we will show 
later, only the monopole receives contributions from all 
super-Hubble modes. 

The power spectrum P^ and amplitude of non- 
Gaussianity, /nl, appearing on the right-hand side of 
Eq.ll]) are those defined in the large volume. How¬ 
ever, looking ahead to the result for the dipole modu¬ 
lation from the local ansatz, Eq.(15l, P^ and /nl will 
be shifted to the observed values. In particular, the lo¬ 
cal non-Gaussianity that generates the observed power 
asymmetry is only the portion that violates Maldacena’s 
consistency relation |54) and is zero in single-clock infla¬ 
tion [^. Since the observed values are ultimately the 


relevant quantities in the analysis, we will not increase 
the complexity of the notation to distinguish the large 
and small volume parameters, except in the appendixes. 

The field (j)(ki) appearing in the integral is no longer 
stochastic but consists of the particular realization of the 
field that makes up the background of a particular Hub¬ 
ble volume. The power spectrum in Eq. © can depend 
on position x within Vs because an individual realiza¬ 
tion (local value) of the fluctuations (/(k^) can be non¬ 
zero. If we consider the average statistics in the large 
volume (equivalent to averaging over all regions of size 
I/s), then the term proportional to /nl in Eq.([^ above 
averages to zero since ((/)(k^))vy = 0. In that case we re¬ 
cover the isotropic power spectrum of Vl- Finally, keep 
in mind that Eq. (§ is still more general than our actual 
GMB sky: it provides the statistics from which to draw 
realizations of our observed modes. Any single sky real¬ 
ization will still be subject to the usual cosmic variance 
that affects the values of small £ modes and that can 
generate a power asymmetry even for volumes where the 
term proportional to /nl is zero. 

It is also interesting to consider a two-field extension 
of Eq.0: 


$(x) = (p(x) -H cr(x) -h /nl^ (o-(x)^ - (o-(x)^)) , (3) 


where both (^(x) and cr(x) are Gaussian random fields 
and are uncorrelated. In this case, the power spectrum 
observed in our sky, P$_s(fc,x), is 


P<i>,s{k,x) = P^{k) 


1 + 4//nLo. 


d?kt 

(27r)3 


a(k,)e**'^'" 


(4) 


where P^{k) = Pipik!) -|- Pa{k) is the mean power spec¬ 
trum in the large volume and / = is the frac¬ 

tion of power in the a field. We will only consider 
cases with weak non-Gaussianity (/nLct^o- ^ 1 and 
Pii>,a ~ Pa)- The amplitude of the local-type bispectrum 
for this weakly non-Gaussian two-field model is given 
by /nl = /^/nLct- Therefore, the inhomogeneous power 
spectrum in terms of the observed /nl and the fraction 
of power / is 


E<j>,s(fc,x) — P^{k) 


, ^/nl [ 


zk^-x 


(5) 

A scale-dependent power fraction /(fc) is a natural way 
to generate a scale-dependent power asymmetry, as can 
be seen from the inhomogeneous term in Eq. 
decreases for large values of k, the modulation of the 
power spectrum will decrease as well. Such mixed- 
perturbation models also have other potentially observ¬ 
able consequences: ^ affects the tensor-to-scalar ratio and 
contributes to large-scale stochasticity in the power spec¬ 
tra of galaxies |55] . 












4 


B. Effect on the CMB sky 

The imprint of the inhomogeneous power spectrum 
given by Eq. (© on the CMB can be described in terms 
of a multipole expansion: 


where jT is a spherical Bessel function of the first kind, 
and X = xn specihes the position of the observed fluctu¬ 
ation: for the CMB, x = rcmb is the comoving distance 
to the last scattering surface. Eq. then implies that 


P^{k,n) = P^{k) 


1 + /nl ^ gLMYLMin ) 


LM 


( 6 ) 


gLM = iSttc 


' |k^|<7r/tc 


(27r)^ 


jL{kix)^{ki)YlM{h)- ( 8 ) 


where is a spherical harmonic, and n is the direction 
of observation on the last scattering surface. To find the 
expansion coefficients gLM we make use of the plane wave 
expansion 

g*k,.x ^ ATTY,i'^jLikex)YlM{k,)YLM{n), (7) 

LM 


The quantity gLM has a fixed value in any single vol¬ 
ume Vs, but when averaged over all small volumes in 
Vl, {gLM)vL — 0 - The expected covariance, on the other 
hand, is non-zero: 


{gLMg*L>M')v^ = 2567r2(-l)^'f^+-^' 
32 


d^kf 


= — Slu^mm' 


= G^TrSLL'dMM' 


l\ke\<-!r/x (^Tr) 

dke k'jj‘l{k(,x)P^{kt)\ 


jL{kix)jL'{kix)P^{kt)Ylj^{ki)YL'M'{ki)\ 


T/a 


dki ,2 


ki 


jL{kix)P^{ki), 


(9) 


( 10 ) 


where in the last line, we have defined the dimensionless 
power spectrum as V^{k) = k^P^(k)/. We have 
again used the subscript Vl to indicate the ensemble av¬ 
erage is over the values of gLM in the full volume Vl- 
Note that both the individual values of gLM and their 
variance depend on the size of the small volume through 
the upper limit of integration in Eqs. Q and While 
the mean statistics in the large volume cannot depend 
on the scale for the small volume, the variance of the 
statistics observed in sub-volumes generically does. It is 
now straightforward to study the monopole and dipole 
contributions from non-Gaussian cosmic variance to the 
modulated component of the power spectrum in a small 
volume. 

In the case of the two-held extension, using Eq. (§ in 
the definition of the modulation moments Eq. gives 




64:TrSLL'dMM' dke .2 


ke 


- -^\9LMgL'M' ^ 


jL{kiX)'Pa{ki)\ 

( 11 ) 


That is, for the same amplitude of non-Gaussianity ob¬ 
served in the 4) held, the variance of the non-Gaussian 
modulations increases by a factor of Icompared to the 
single source = I) local model. 


1. Monopole modulation (L = 0) 


The power-spectrum amplitude shift, Aq, in the 
parametrization of Eq. is: 


P^{k)=P^{k) [I + Ao] 


pm 


1 + /nl 



( 12 ) 


where Aq can be either positive or negative, but has a 
lower bound Aq > —I. From the discussion above and 
Eq.(|^, it is clear that goo is Gaussian distributed with 
zero mean and variance given by: 


(goo) = 647r 


f dke 

sirL{kex) 

J ke 

kex 


P<j>{ke). 


(13) 


Therefore, the distribution of the monopole power mod¬ 
ulation amplitude Aq also follows a normal distribution, 
for small values of Aq, (|^o| ^ 1)) with zero mean and 
standard deviation: 

(14) 


The expression for (gpo) i® sensitive to the infrared limit 
of the integral. That is, all super-Hubble modes can con¬ 
tribute. Interesting aspects of cosmic variance arising 
from this term, including effects on the observed non- 
Gaussianity in small volumes have been subjects of in¬ 
vestigation in [M1162) . In particular, the observed value 
of /nl is, in general, shifted from the mean value in the 
large volume. 
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For a constant /nl, the effect of the monopole modu¬ 
lation is to change the power-spectrum amplitude on all 
scales and therefore is not observationally distinguish¬ 
able from the “bare” value of the power-spectrum am¬ 
plitude. For scale-dependent non-Gaussianity, there is 
a scale-dependent power modulation, which can generi- 
cally be interpreted as shifting the spectral index in the 
small volume away from the mean value in the large vol¬ 
ume. In cases where the amplitude of non-Gaussianity 
is small (and consistent with zero) at small scales (large 
£), the power-spectrum amplitude from those scales can 
be taken as P^(k), and then one can look for monopole 
modulation at large scales for which the non-Gaussianity 
constraints are not as strong. The large-scale power sup¬ 
pression anomaly [531 El] is exactly such a situation. We 


will return to this point in more detail in Section IV 



modulation multipole, L 


2. Dipole modulation (L = 1) 

The dipole modulation of the power spectrum in the 
parametrization of Eq. <§ is given by: 


P$(fc,h) = P4,{k) 


1 +/nl ^ giMYiM { n ) 
M=- 1 , 0,1 


(15) 

Since we are interested in the dipole modulation of the 
observed power spectrum in the CMB sky, the above 
equation should be obtained from Eq.([^ by absorbing 
the (unobservable) monopole shift to the observed power 
spectrum. Then, on the right-hand side of Eq.([T5|), /^(A:) 
is the observed isotropic power spectrum and jnl is the 
observed amplitude of local non-Gausianity within our 
Hubble volume. See Eq. ( |B11 ) and the discussion there 
for details. (Appendix [Cj contains the corresponding ex¬ 
pression in terms of bipolar spherical harmonics.) The 
giM coefficients are Gaussian distributed with zero mean 
and a variance 


{giHglm) = 647r 


f dke 

sin(fc^a:) 

cos{kix) 

I h 

{ktxY 

kiX 




If we pick a direction di in which to measure the dipole 
modulation Ai such that 


P<s,{k) = P^{k) [1 -I- 2Ai COS0] 


(16) 


where cos 9 = di ■ n, then the contribution to the dipole 
from the non-Gaussianity is = jy^/NL^io, which 
is normally distributed with mean zero and standard de¬ 
viation: 


O'/nl = W^I/NLKffm)^ 


(17) 


In the two-field model Eq.(|^, using Eq.([II|), the standard 
deviation gets modified: 


'^/nl ['^/ni.]{=1 • 


(18) 


FIG. 1. (color online). The expected higher order modula¬ 
tions for /nl = 500,250 (red square, open green square). We 
see a sharp decrease in the expected modulation amplitude 
from superhorizon fluctuations at higher L multipoles. The 
blue circle for L = 1 is the variance measured from our Gaus¬ 
sian Sachs-Wolfe CMB maps. The dependence of the corre¬ 
sponding variances for higher multipoles L in Gaussian CMB 
maps on the multipole L is much weaker than the sharp falloff 
seen in the expected amplitude due to local non-Gaussianity 
(See Figure 2(d) of [ 6 ], for example). 


This shows that for ^ < 1 (e.g. a mixed inflaton-curvaton 
model), it is easier to generate the hemispherical power 
asymmetry with a small value of /nl- However, there is 
a minimal value of ^ that can generate a power asymme¬ 
try of a given amplitude: the requirement of weak non- 
Gaussianity in the non-Gaussian field [cr(x) in Eq. @] 
demands that ^ 


The above discussion of the distribution of the dipole 
asymmetry A^ assumes that we measure A^ in a fixed 
direction di. However, we have no a priori choice of di¬ 
rection di in most situations. This is especially true when 
considering a power asymmetry that is generated by the 
random realization of superhorizon perturbations as op¬ 
posed to a single exotic perturbation mode. Therefore, 
observations of dipole power modulations are necessar¬ 
ily reported using the amplitude of dipole modulation 
in the direction of the maximum modulation. To ob¬ 
tain that amplitude, we can consider any three orthonor¬ 
mal directions {di,d 2 ,d 3 ) on the CMB sky and measure 
the corresponding three dipole modulation amplitudes 
(Ai, A 2 , A 3 ) in the three corresponding orthonormal di¬ 
rections for each sky. The amplitude of modulation for 
the CMB sky (simulated or observed) is then given by 
A = {A\ -I- A 2 -I- A |)5 . Clearly then, A follows the y dis¬ 
tribution with three degrees of freedom (also known as 
the Maxwell distribution). In the Section HI A we will 
directly test the distributions and parameters obtained in 
this section using numerical realizations of CMB maps. 
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3. Higher multipole modulations 


The anisotropic modulation of the power spectrum is 
expected to continue to higher multipoles in the presence 
of non-Gaussianity. However, as shown in Figure the 
expected value of the modulation gets smaller quickly for 
higher multipoles L. The corresponding expected vari¬ 
ance of higher multipole modulations for Gaussian CMB 
maps, however, is only weakly dependent on L. See, for 
example. Figure 2(d) of [S]. There is no evidence for mod¬ 
ulation at higher order multipoles in the Planck tempera¬ 
ture anisotropies data (see Figure 34 of HI)- In statistical 
analysis of the kind we discuss later in Section |IIID[ it 
may, nevertheless, be useful to add higher multipole mod¬ 
ulations (at least the quadrupole L = 2) at large scales 
as it may provide increased evidence for or against non- 
Gaussian mode coupling. An approximate constraint on 
/nl 32 M may be obtained from the result for the Fourier 
space quadrupole modulation constraint in [65] • There 
are two possible scenarios: (i) the expected amplitude 
of the modulation is larger than that from cosmic vari¬ 
ance in the Gaussian case; in this scenario, the lack of 
observation of such a modulation in the data will disfa¬ 
vor the non-Gaussian model that is used to explain the 
monopole and dipole modulations, (ii) The expected am¬ 
plitude of the modulation is within the cosmic variance 
from the Gaussian case, in which case the data are not 
discriminatory for or against the non-Gaussian model. 


C. Connection to prior work 


Before proceeding, we pause to connect Eq.(15) to the 


EKG mechanism to better see how scale-dependent non- 
Gaussianity eliminates the need for enhanced perturba¬ 
tions on superhorizon scales. In the EKC mechanism, 
a single superhorizon perturbation mode in a specta¬ 
tor held during inflation is responsible for generating 
the asymmetry; the original proposal used a curvaton 
held m, but later work extended the mechanism to any 
source of non-Gaussian curvature huctuations [HI l34] . 
For example, consider a held a that generates a curvature 
perturbation C. = JV^Sa. A superhorizon (SH) sinusoidal 
huctuation in tr. 


ctsh(^) = ctl cos ■ x + o'j 


(19) 


will generate a dipolar power asymmetry in the curvature 
power spectrum [34]: 


Pcik,x) = Pi;{k) 


12 

1- ^hhikL ■ x)NaCrL sin 9 

5 


( 20 ) 


to hrst order in k^x. In terms of the Bardeen potential, 
4> = (3/5)^, the power asymmetry is 


P,i,{k,x) = P^{k) 1-4/nl(A:l •5')4‘Lsin6> , (21) 


where = {'i/'5i)NaaL- The Fourier transform of the 
superhorizon huctuation given by Eq. (19) is 


®Sh(^ = (27r) 




P doikh - ^ + e 5D(kL+^ 


( 22 ) 

Inserting this expression into Eq.® for qlm implies that 


gw = -167r$^ 


k^xsinO -\- 0{kj^x^). 


(23) 


With this expression for gio, Eq. (15) matches Eq. (21) 


to hrst order in k^x. Therefore, we see that the EKG 
mechanism can be described by our framework. 

For a single superhorizon mode, Eq. ( |21| ) implies that 
the non-Gaussian contribution to the dipole is 


= 2|/nl|A$, 


(24) 


where A<I> = ^Lk^xsinO is the variation of $ across the 
surface of last scatter. Since k^x < 1, > A$. The 

RMS amplitude of 4) values given by extrapolating the 
observed value of Vc^ to larger scales is 




— \/^) 


(25) 


where 7^,^ ~ 8 x 10 [55]. It follows that the amplitude 

of the superhorizon mode is bounded from below as 


^rn 


> 


A 


2 |/nlIx/^' 


(26) 


If A = 0.06 (and is entirely due to the non-Gaussianity) 
and |/nl| < 100, then > 10.6, which implies 

that the superhorizon mode must be at least a lOcr huc¬ 
tuation. This is why Eq.( [19[ ) was not originally consid¬ 
ered to be part of the inhationary power spectrum, but 
rather a remnant of pre-inhationary inhomogeneity or 
a domain-wall-like feature in the curvaton held m- It 
was then necessary to consider the imprint this enhanced 
superhorizon mode would leave on large-scale tempera¬ 
ture anisotropies in the GMB through the Grishchuk- 
Zel’dovich (GZ) effect [57]. Although the curvature per¬ 
turbation generated by Eq.(19) does not generate an ob¬ 


servable dipolar anisotropy in the GMB [SMO], it does 
contribute to the quadrupole and octupole moments, and 
observations of these multipoles severely constrain mod¬ 
els that employ the EKC mechanism mm- 

However, if we relax our upper bound on /nl to 270 
or 500, Eq.(26) indicates that a Aa or 2 ct huctuation, re¬ 


spectively, could generate an asymmetry with A = 0.06. 
The odds of generating the observed asymmetry are also 
improved by accounting for the fact that there are three 
spatial dimensions, which provide three independent op¬ 
portunities for a large-amplitude huctuation. We will 
see in Section [HI Dj that considering the combined con¬ 
tributions of several superhorizon modes and accounting 
for the red tilt of the primordial power spectrum fur¬ 
ther increases the probability of generating the observed 
asymmetry, to the point that the p—value for A = 0.06 
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increases to greater than 0.05 for |/nl| ^ 300. Thus, 
if the perturbations on large scales are sufficiently non- 
Gaussian, there is no need to invoke enhanced superhori¬ 
zon perturbations to generate the observed power asym¬ 
metry. 

In the absence of an enhancement of the superhorizon 
power spectrum, the variance of the quadrupole moments 
and octupole moments in the CMB will not be altered. 
Consequently, we do not expect significant constraints on 
such models from the GZ effect. We note though that the 
specific realization of modes outside our sub-volume will 
still source quadrupole and octupole anisotropies in the 
CMB. For realizations that generate a large power asym¬ 
metry, the GZ contribution to these anisotropies would 
likely be larger than expected from theoretical predic¬ 
tions of C 2 and C 3 and aligned with the power asymme¬ 
try. However, this effect may be difficult to disentangle 
from the monopole power modulation described in Sec¬ 


tion IIB 1 and we leave a detailed analysis of this obser¬ 
vational signature to future work. 

Furthermore, the ratio for a given value of 

A and /nl can be significantly reduced if we consider 
mixed Gaussian and non-Gaussian perturbations. Using 
the mixed perturbation scenario introduced in Eq.([^, 
(with ^ 


$ 


L.u 


A 

^ r 


> 




> 


^3/2 


2|/nl|v^ 2|/nl|v^' 


(27) 


In the last inequality, we employ the fact that ^ > A 
is required for the non-Gaussianity in the (t(x) field to 
be weak enough to make the ©(/nl^T’o-) contribution to 
P^^a{k) negligible. In this case, could be sourced by 
a Icr fluctuation in the cr(x) field if |/nl| = 270. The pos¬ 
sibility of using a mixed curvaton-inflaton model to gen¬ 
erate a scale-dependent asymmetry using a single large 
superhorizon perturbation was explored in |12j . 


III. STATISTICAL ANISOTROPY IN THE CMB 
POWER SPECTRUM 

A. Numerical tests 


We now present numerical tests of our analytic expres¬ 
sions for the dipole power modulation in the case of local 
non-Gaussianity. We will work in the Sachs-Wolfe (SW) 
regime: we only consider 


AT _ $ 


(28) 


We generate 10000 simulated Gaussian SW CMB skies 
using 


Att dh 

= (30) 

for 0 < £ < 300; the primordial power spectrum is given 
by 


( , \ ria —1 

yJ (31) 

with A^ = 7.94 x 10“^° and Ug = 0.965 (from Planck 
TT,TE,TE-|-lowP column in Table 3. of [ 66 ]), and ko = 
0.05 Mpc“^ as the pivot scale. 

Then it is easy to generate non-Gaussian Sachs-Wolfe 
CMB temperature maps using Eq.(|2^ for a constant 
/nl- Unlike most CMB analyses, we will keep the dipole 
variance term Ci. A non-zero Ci is used to model the 
dipolar anisotropy in density fluctuations on the scale 
of the observable universe (from the perspective of the 
large volume V/). However, note that the Ci we use 
is not what we would measure for the CMB dipole, 
even if we assume that the dominant contribution to 
the measurement of the dipole from our local motion 
m has been subtracted out. This is because, for adi¬ 
abatic fluctuations, the leading-order contribution to the 
observed CMB dipole from superhorizon perturbations 
exactly cancels the Doppler dipole generated by the su¬ 
perhorizon perturbations [69l|70|. 

It is convenient to set the monopole Cq to zero for the 
purpose of studying dipole modulations; otherwise, the 
cosmic variance power asymmetry (i.e. the contribution 
that is not due to local non-Gaussianity) will be differ¬ 
ent for the weakly non-Gaussian realization compared 
to the Gaussian realization from which it is generated. 
Therefore, in this section, we use numerical realizations 
with non-zero Cq values only when testing the monopole 
modulation formula. The expression for Cq is infrared di¬ 
vergent, so we assume an infrared cutoff fcmin; the same 
cutoff scale is used to compute the expected amount of 
monopole power modulation Aq. Numerically, 



sin kx 
kx 


(32) 


The fcmin cutoff Can be related to the number of super¬ 
horizon e-folds of inflation (if interpreted as such) as: 
-/Vextra = In [(7r/rcmb)/^min] ■ The above integral gets 
most of its contribution from k < Tr/rcmb) and therefore 
can be well approximated by jSSj: 


Therefore, for local non-Gaussianity, Eq. the temper¬ 
ature fluctuation is given by: 


AT 


/nl 


AT 

~Y 


gaus 


3/nl 


AT 

~Y~ 


2 

gaus 


AT 

~Y~ 



Cq 



/ \ ris —1 

/ ^ \ 

'1 _ g-(n.-l)Arextra' 

V ^O^cmb J 

rig - I 


(33) 


for Us yf 1, and ( 47 r/ 9 )A 0 Aextra for Ug = 1. Additional 
details about our numerical results can be found in Ap¬ 
pendix [^ 



















B. Monopole modulation (L = 0) 


The normally distributed monopole shift amplitude Aq 
for a local non-Gaussian model is given by Eq.(14), and 
can be written in terms of Co using Eq.(32| as: 


^/nl 


6|/nl| 


(34) 


The necessary infrared cutoff has already been set by 
the value of fVextra in Eq.(33) to compute Cq. Although 


it is not possible to observe monopole modulations for a 
constant local /nl, we can test the expected modulations 
assuming a value of A^extra- The probability distribution 
of the shift for any /nl 7 ^ 0 is 


PN{Ao,a^°n 


1 

cr“™°v^ 


Exp 


1 ( -^0 

2 (^(jinono 


(35) 


where the variance has contributions from the Gaussian 
realization and the non-Gaussian coupling to the realiza¬ 
tion of long wavelength modes: ((j™™°)2 = (cr^°"°)^ -I- 
(crQ™°)^. For our numerical tests, is the variance 

of Aq measured in Gaussian GMB maps. The quantity 
we measure for Aq from each realization of CMB maps 
is: 


An = 




£( 2 ^+ 1 ) 


Ce - Cf 

^true 


(36) 


where are the input angular power-spectrum values 
used to obtain the set of numerical CMB maps, and Cg 
is the angular power spectrum of a particular realization 
of that set of GMB maps, and £max = 100. In Figure 
we plot the distribution of Aq, using Eq.([^, for /nl = 
0, 50,100, along with the distribution obtained from the 
numerically generated Sachs-Wolfe CMB maps. 


C. Dipole modulation {L = 1) 


-^extra=50 



Ao 


FIG. 2. (c olor online). Test of the monopole modulation 
formula Eq.( |35[ ) for the local non-Gaussian model with /nl 
specified in the figure. The dotted blue line for the /nl = 0 
(Gaussian) model is the best-fit normal distribution to the 
distribution obtained from numerically generated Sachs-Wolfe 
GMB maps, while the other two curves (dashed green, /nl = 
50 and solid red, /nl = 100) are obtained using Eq.(35l. 
The value of the CMB monopole Co is set using Eq. (331 with 
A^extra = 50. The measurement of Aq is bounded below i.e. 
Ao > —1. The normal distribution is an excellent fit for 
small modulations i.e. |Ao| <C 1, whereas the distribution is 
positively skewed for larger values of |Ao|. 


non-Gaussian maps, the distribution of the A^s have an 
increased variance given by: cr^ = + CTq, where 

is given by Eq.( |17[ ). The power asymmetry dipole am¬ 
plitude for each CMB sky is then A = (A^ + A 2 + A|) ^. 
The probability distribution function (pdf) of A is the x 
distribution (or the Maxwell distribution): 


Px(^-cr) 



■-A 2 - 

2 (t 2 


(38) 


A dipole modulation of the power spectrum defined as 
in Eq. (16 1 generates a hemispherical power asymmetry 
with the same amplitude 2Ai. Therefore, we will look at 
the quantity: 


A, = 


Etr(2^ + i) 


£( 2 £ + 1 ) 


e =2 


AC, 

2C, ’ 


(37) 


with /niax = 100 and AC, = C/ — C/, where -I- and — 
refer to two hemispheres in some direction di. We will 
consider AjS in three orthonormal directions di, d 2 , and 
da on the sky. Each Ai in a particular direction di is 
normally distributed with zero mean. The variance 
can be measured from the numerical realizations of Gaus¬ 
sian Sachs-Wolfe CMB maps and depends on the CMB 


multipoles used in Eq.(37l and the value of the C,s. For 


shows that the distri- 
jtained from the CMB 


where a = -|- a^. Figure 

bution of asymmetry amplitudes o' 
realizations agree extremely well with the x distribution 
given above. Note that only ao is measured from the 
numerical maps; is directly computed for a value of 
/nl using Eq. Q- 


D. Statistical analysis 

In this section, we present examples of how we can per¬ 
form a statistical analysis using the results from the pre¬ 
vious section for the distributions of power modulation 
on the Sachs-Wolfe CMB sky. While direct comparison 
of the amplitudes obtained in our Sachs-Wolfe CMB re¬ 
alizations with the reported values of power asymmetry 
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FIG. 3. (color online). Left: The distribution of power asymmetry Ai (in a particular direction di) measured in 10000 simulated 
CMB skies as described in the text. From each simulated map, three Ai values are generated in three orthonormal directions in 
the sky. The dotted blue line for the Gaussian GMB maps is the best-fit normal distribution curve; this gives us the Gaussian 
cosmic variance standard deviation ao- The curves (dashed green and solid red) for the non-Gaussian mode ls are normal 
distributions with zero mean and variance given by = ctq -|- where is computed using Eq.(17|. Right: The 

distribution of the amplitude of power asymmetry A from the simulated Gaussian and non-Gaussian GMB maps. The lines are 
the computed Maxwell distributions for corresponding o values and match the distributions obtained in numerical realizations 
very well. Note that the figures above were generated for the single source local model. If we used the mixed infiaton-curvaton 
two-field extension with the curvaton power fraction ^ < 1, we would get the above distributions for smaller values of /nl- For 
example, for ^ = 0.25, the distributions shown for /nl = 500 above would be generated by a smaller /nl = VO.25 x 500 = 250. 


A is not possible, we can make a connection between our 
simpler case and the asymmetry in the observed CMB 
sky by using the p-value of the asymmetry. For a given 
measurement of A and the normalized pdf for A, p{A) 
[which in our model depends on /nl, see Eq.(38)], the 
p-value is simply given by p{A')dA'; i.e. it gives the 
probability that the observed value of the asymmetry am¬ 
plitude is greater than some threshold value A. We find 
that an asymmetry amplitude oi A = 0.055 is approxi¬ 
mately 3.3 ct, i.e. a p-value of 0.001 with respect to the 
distribution of A obtained in our Gaussian Sachs-Wolfe 
CMB maps. This is approximately equal to some of the 
more recent reports for the significance of the hemispher¬ 
ical power asymmetry 0 [72]. Therefore, we will use 
A — 0.055 as the value of the asymmetry when making 
connections with the observations of the anomaly. 

When we have a measurement of the power asymme¬ 
try amplitude A, we can write the likelihood for /nl 
as £(/nl|^) = O’), whose expression is given in 

Eq.([M|. From this likelihood, we can infer the posterior 
distribution for /nl given a measurement of A. We can 
interpret the statistics in different ways: 


• We can use any power asymmetry as a signal of 
local non-Gaussianity. Using only the large-scale 
CMB multipoles {I < 100), for a given value of A, 
we can obtain the posterior distribution for |/nl| 
(averaged) for the corresponding range of scales. In 
Figure]^ we plot the /nl posterior for a few values 
of the asymmetry A, assuming a uniform prior on 


I/nl|- 

• We can combine the large-scale bispectrum con¬ 
straints on /nl with the constraints from the power 
asymmetry A. For this, we use a rough estimate 
of /nl for £ < 100 of /nl = -100 ± lOO(lcr) (es¬ 
timated from Figure 2 of [S2|). We assume that 
the /nl posterior from WMAP is a normal distri¬ 
bution. However, since the power asymmetry is 
only sensitive to the magnitude of /nl and not the 
sign, we use the folded normal distribution (given 

/measured _ ,, _1_ 

/nl — /X It aj. 


7\/2‘K 


Pfold(/NL) = < 


. 1 ( /nl-^ 


. 1 / X’nl+m 


/nl > 0 
/nl < 0. 


Then, we multiply the above pdf with £(/nl|A) to 
get the combined likelihood from which we can get 
the posterior for /nl after normalizing. We show 
the combined posterior distribution of /nl for a few 
power asymmetry amplitudes A in Figure]^ 

• Although the power asymmetry data alone {A « 
0.055) prefers |/nl| ~ 500 as the most likely value 
(or |/nl| ~ 200 when the bispectrum constraints 
are also applied), the probability of an asymme¬ 
try increases whenever /nl 0. In Figure we 
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FIG. 4. (color online). The posterior probability distribu¬ 
tion of |/nl| values for different observed amplitudes A of 
power asymmetry. Only large-scale CMB modes I < 100 are 
used to compute the cosmic variance pdf for A. A = 0.055 
corresponds to a p-value of 0.001 in our /nl = 0 numeri¬ 
cal maps. This is about 3.3 (j, approximately equal to some 
of the reported significance of the power asymmetry anomaly 
[ 61122 ]. Although our formula for the expected asymmetry be¬ 
comes less accurate for larger values of /nl we have checked 
that it approximates the numerical results quite well even for 
/nl = 2000. Therefore, the shape of the posterior distribu¬ 
tions obtained above in the /nl window shown will not be 
affected by the inaccuracy of the formula at larger /nl val¬ 
ues. However, the change in the distribution for larger /nl 
can change the normalization. 


quantify how the probability of an observed dipole 
modulation changes as |/nl| increases. For exam¬ 
ple, the p-value of 0.001 for A = 0.055 changes to 
0.046 for |/nl| = 265 (which is within the 2cr win¬ 
dow of the large-scale bispectrum constraint). In 
other words, even an amplitude of non-Gaussianity 
well below |/nl| ~ 500 renders the observed asym¬ 
metry less “anomalous.” 


E. Bayesian evidence 


The previous section demonstrated that the power 
asymmetry data can be used to constrain non-Gaussian 
models, and that the amplitude of the observed asym¬ 
metry is less “anomalous” when non-Gaussianity is in¬ 
cluded. However, we also need to ask whether the data 
are such that the non-Gaussian model is preferred over 
the Gaussian. 

To compare the posterior odds for different models A4i, 
given the data y, we compute the Bayes factor 


pjylMi) 

piy\M2)’ 


(39) 


FIG. 5. (color online). The posterior probability distribution 
of |/nl| after combining the bispectrum constraints at large 
scales (Z < 100, /nl = —lOOTlOO) with the power asymmetry 
constraint for the given value of A. 


where the factors in the numerator and denominator are 
the model likelihoods for models 1 and 2 respectively. In 
the simplest comparison, we take /nl as the only pa¬ 
rameter of the models. The data we consider include 
the measured amplitude of the power asymmetry and the 
CMB constraint on the amplitude of the local bispectrum 
on large angular scales. For an introduction to Bayesian 
statistical methods applied to cosmology, see for example 

m- 

The non-Gaussian model reduces to the isotropic 
Gaussian model for /nl = 0 (while the probability of 
the power asymmetry remains non-zero). In that case 
the evaluation of the Bayes factor can be simplified and 
a direct Bayesian model comparison can be done using 
the Savage-Dickey density ratio (SDDR) [731 IZl]- The 
SDDR is given by. 


p{&z\y,Mi) 

7r(6»i|A4i) 


(40) 


where A4i is the more complex model (non-Gaussian 
in our case) that reduces to the simpler model A4o 
(Gaussian) when the set of parameters 9i goes to 9* 
(/nl -f /nl = 0). Here, p{9i\y,Mi) is the posterior for 
/nl (plotted in Figure and 7r(0i|A4i) represents the 
prior for the parameter in the complex model A4i. Our 
current case only has one parameter (/nl) and one datum 
(the dipolar asymmetry A) . While there may be other in¬ 
teresting possibilities to consider for the prior probability 
of /nlj we illustrate the calculation of Boi above using 
the constraint on the parameter /nl from large-scale bis¬ 
pectrum measurements as reported by the WMAP and 
Planck missions as the prior. 

In Table |Tj we list SDDR for a few values of Aobs- For 
the prior, we have used the folded normal distribution for 
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FIG. 6. (color online). The p-value for different values of 
asymmetry amplitudes A (i.e. the probability of obtaining an 
asymmetry amplitude equal to or greater than A) in a local 
non-Gaussian model as a function of the value of |/nl|. For 
A = 0.055, we see that the signihcance goes below 3cr around 
|/nl| ~ 100. For the two-field extension of the local model 
(with a curvaton power fraction described in the text, the 
x-axis should be labeled |/nl|/v^- Then, for ^ < 1, a smaller 
bispectrum amplitude |/nl| (by a factor of than labeled 
in the figure above) is required to achieve a p—value shown in 
the figure above. 

/nl = —100 ± 100 which is a rough estimate of /nl for 
the largest scale i.e. up to £ = 100 from (52]. Note that 
the only value from the prior pdf that is used to compute 
the SDDR is /nl = 0, so the above consideration from 
large-scale /nl constraints is the same as using a uniform 
prior for |/nl| in the range (|o, ~ (0,207). If 

the prior range is expanded, then the magnitude of i?oi 
increases thereby reducing the evidence for non-zero /nl- 
For ^obs = 0.055 (whose p-value roughly corresponds to 
observed A ), the strength of evidence for a non-zero /nl 
is between weak and moderate in the empirical (Jefferys’) 
scale m quoted, for example, in Table 1 of m- 

The results in Table I show that the data we have 
used, at least in this simple analysis, show no more than 
a weak preference for the non-Gaussian model. A more 
thorough analysis is unlikely to change this conclusion 
very much: in m, the authors use earlier studies of 
the power asymmetries in the WMAP data to put the 
best possible Bayesian evidence of InBpi « —2.16 corre¬ 
sponding to odds (< 9 : 1, weak support). The method 
to compute the maximum possible Bayesian evidence is 
based on Bayesian calibrated p-values m\- The p—value 
used from the data analysis of the 3-year WMAP maps 
was p = 0.01 m- If one instead used p = 0.001, which 
is approximately the level of significance from various 
more recent analyses of Planck and WMAP tempera¬ 
ture anisotropy maps, the best possible Bayesian evi¬ 
dence in favor of the A ^ 0 anisotropic model becomes 


-^obs 

p-value 

SDDR (Boi) 

InSoi 

0.02 

0.5511 

1.1362 

0.128 

0.04 

0.0381 

0.6174 

-0.482 

0.05 

0.0043 

0.3211 

-1.136 

0.055 

0.001 

0.2012 

-1.603 

0.06 

0.0003 

0.1123 

-2.186 


TABLE I. Savage-Dickey Density Ratio (SDDR) for different 
observed values of dipole power modulations at large scales. 
The p—values listed above are computed with respect to the 
distribution of A values from Gaussian GMB maps. A reason¬ 
able value to compare to various reports of the observed hemi¬ 
spherical power asymmetry is a p—value of 0.001 (« 3.3a) i.e. 
A = 0.055. Values of |lnBoi| = 5.0, 2.5 suggest strong and 
moderate evidence, respectively m- In our convention above, 
a negative value for the logarithm of the Bayes factor means 
the evidence is in favor of the more complex non-Gaussian 
model All. 


lni?oi ~ —4.0 corresponding to the odds (< 50 : 1). A 
p—value of 0.0003 (about 3.6cr) is necessary to obtain a 
best possible InRoi ~ —5.0, which implies strong evi¬ 
dence [7^ . 


IV. BEYOND THE LOCAL ANSATZ 

The previous section discussed in detail the effect of 
local-type non-Gaussianity with constant /nl that cou¬ 
ples a gradient (induced by superhorizon modes) across 
the GMB sky to the observable modes. We demon¬ 
strated that a dipolar asymmetry is expected in models 
with local-type non-Gaussianity. Of course, local non- 
Gaussianity as the source of the asymmetry is only com¬ 
patible with the data if we restrict ourselves to the largest 
scales. Both the amplitude of the asymmetry and the 
amplitude of non-Gaussianity must sharply decrease on 
smaller scales. The question then is whether there is a 
different model of non-Gaussianity that is consistent with 
all observational constraints and generates the observed 
asymmetry in detail. If so, does the current data fa¬ 
vor this model over the isotropic, Gaussian assumption? 
Gould future data ever favor such a model? 

To address some of these questions, we will first demon¬ 
strate that scale-dependent modulations are a generic 
feature of non-Gaussian models other than the local 
model. We will then construct a scenario that is more 
likely to be preferred by the data by considering scale- 
dependent local non-Gaussianity. Finally, we will pro¬ 
vide examples of evaluating the Bayesian evidence for 
this scenario. Although a model of non-Gaussianity be¬ 
yond the local ansatz may add more parameters, if the 
model has other consequences in the data we might hope 
to find more evidence for it. This is particularly true if 
the model has measurable effects on smaller scales, where 
the usual cosmic variance for Gaussian models is smaller. 
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A. Power asymmetry from general bispectra 


We can easily extend the inhomogeneous power spec¬ 
trum calculation in the presence of local non-Gaussianity 
to other bispectrum shapes. For example, consider that 
a Fourier mode of the Bardeen potential is given by |62j : 

$(k) = (j){k) +'^ J J 0(qi)</'(q2) 

^2(qi,q2,k)i5D(k - qi - q 2 ) + • ■ -(41) 

where as before (j) is a Gaussian field. The kernel N 2 
can be chosen to generate any desired bispectrum and 
the dots represent terms higher order in powers of (j) 
(which generate tree-level n-point correlations). Gonsid- 
ering only the generic quadratic term, the power spec¬ 
trum in sub-volumes can be computed as in the case 
of the local bispectrum (see Appendix and we get 
P$,s(fc,x): 




1 + 2/nl 


(27r)3 


<).(k,)iV2(k,,-k,k)P'^^’' 


(42) 

From the form of the above equation, one can see that a 
fc-dependent power modulation is a feature of non-local 
non-Gaussianity i.e. the k dependence of the kernel N 2 is 
carried by the modulated component of the power spec¬ 
trum in the small volume. The kernels for local, equilat¬ 
eral and orthogonal bispectrum templates are [62]: 


^local ^ 2; 

Arortho 4A;| - 2kkg _ 

,requil _ 

^^2 - ^2 • 


(43) 


If one uses the kernel for equilateral- or orthogonal- 
type non-Gaussianities, then the monopole shifts are not 
infrared divergent. However, the magnitudes of modula¬ 
tion (both monopole and dipole) are smaller compared to 
the local case i.e. a very large amplitude of or /nl^° 
is necessary for the effect to be interesting. For example, 
we plot the expected modulation amplitude for local- 
, equilateral- and orthogonal-type non-Gaussianities in 
Figure In Figure we illustrate that for local, or¬ 
thogonal and equilateral bispectra, the power asymmetry 
is generated by perturbation modes that lie just outside 
the horizon. The quasi-single held model m may also 
be interesting to consider: it has a scale-independent bis¬ 
pectrum with a kernel that varies between the local and 
equilateral cases depending on the mass of an additional 
scalar held coupled to the inhaton. 

In general then, if the power asymmetry is coming 
from mode coupling, the fact that the observed asym¬ 
metry falls off on small scales implies that shorter scales 
are more weakly coupled to superhorizon modes than 
larger scales are. This is possible with either a scale- 
independent bispectrum (as the equilateral and orthogo¬ 
nal cases above demonstrate) or with a scale-dependent 
bispectrum. 



FIG. 7. (color online). The expected dipolar modulation of 
the power spectrum for large amplitude local-, orthogonal- 
and equilateral-type non-Gaussianities. While the amplitude 
of expected modulation is smaller for non-local shaped bispec¬ 
tra, the modulation generated by them is scale dependent. 



FIG. 8. (color online). The expected amplitude of the power 
asymmetry from superhorizon modes {kt < vr/rcmb), as a 
function of the minimum wavenumber fcr.min considered to 
compute for large-amplitude local-, orthogonal- and 

equilateral-type non-Gaussianities. We can see that most of 
the contribution is due to the modes with wavelengths al¬ 
most the size of the observable universe i.e. within one e- 
fold. We get similar behaviors for the monopole modulations 
(not shown in the figure) except for the case of local non- 
Gaussianity. For local non-Gaussianity, the monopole modu¬ 
lation gets contributions from arbitrarily small k modes and 
(S'oo)iocai becomes infrared divergent. 
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B. Generating a scale-dependent power asymmetry fNL{k): 


To match the observed scale dependence of the power 
asymmetry anomaly, the strength of coupling of subhori¬ 
zon modes to the long wavelength background must be 
scale dependent. The relevant scale dependence in this 
context can be fully parametrized by introducing two bis- 
pectral indices that capture the scale dependence in our 
observable volume and a more general coupling strength 
to the long wavelength modes: 


P‘S>,sik,^) 


pm 


l+4/NL(fco)(^) ' 


X 


(27r)3 





(44) 


Here n/ < 0 turns off any power asymmetries on shorter 
scales. The parameter a < 0 enhances the sensitivity 
of the model to infrared modes (as used in [131ISI])- In 
the case a < — 1, the dipole asymmetry would be in¬ 
frared divergent in a universe with a scale-invariant or 
red-tilt power spectrum. Notice that scale-invariant bis¬ 
pectra always have Uf = —a, so any scale-invariant bis¬ 
pectrum that increases IR sensitivity also increases the 
expected asymmetry on smaller scales. Finally, although 
we have used /nl to label the coefficient above, a simi¬ 
lar expression can be derived from higher-order correla¬ 
tion functions (e.g., to capture the effects of ^nl [S3])- 
Previous discussions of the power asymmetry from scale- 
dependent non-Gaussianity include [1211131 UHl ESI l30l - 

isa- 

In principle, additional data could eventually constrain 
all of the parameters introduced above (or at least their 
values on subhorizon scales). However, here we will con¬ 
sider only one additional measurement (the large-scale 
power suppression) and so we will restrict our attention 
to the case with just one additional parameter. We take a 
local-shape bispectrum with an amplitude that depends 
on the scale of the short wavelength mode as 


/NL(fc) = /^L ( ^ 


(45) 


In terms of the parameters in Eq. (441, n/ = n and 
a = 0. 

In addition to a scale-dependent power asymmetry 
and a scale-dependent bispectrum amplitude, a scale- 
dependent local non-Gaussianity also generates a scale- 
dependent modulation of the power-spectrum amplitude. 
This is the monopole power modulation (L = 0) dis¬ 
cussed in Section |HB (a similar point was made in 


[30j ). When local-type non-Gaussianity has a scale- 
independent amplitude, the power spectrum amplitude 
is modulated similarly at all scales, thereby making the 
effect unobservable. However, the scale-dependent case 
is more interesting as it makes the power modulation 
scale dependent and therefore an observable effect. This 


can be easily seen by generalizing Eq.(I2) for the case of 


P$(fc) = Prj,{k) 


1 + /NL(fc) 


goo 

2^/n 


(46) 


As previously discussed, goo is normally distributed with 
zero mean and the variance (goo) requires a cutoff to 
limit contributions from arbitrarily large modes, which 
we have earlier parametrized as the number of super¬ 
horizon efolds of inflation. While the above formula is 
only valid for small modulations Ao, Figurej^shows that 
the formula is quite accurate at least up to |Ao| « 0.3. 
There is an additional subtlety because, in the presence 
of a monopole shift, the observed bispectrum on large 
scales will not be exactly of the form in Eq.(45|. How¬ 
ever, the difference is small for small Aq. 

Gonsider a simple example of scale-dependent lo¬ 
cal non-Gaussianity given by = —0.64, ko = 

60(7r/rcmb) and /^l = 50. In multipole space, one can 
approximate /^l = /nl(^ = 60) = 50 and /nl(^) = 
50(^/60)“3'3‘^. These numbers are chosen to facilitate 
comparison with the scale-dependent modulation model 
results in [9] (see Table I therein). In Figure]^ we plot 
/nL; (tke expected Icr amplitude of the power asym¬ 
metry in a particular direction due to /nl), and 
(the expected Icr shift in the amplitude of the power spec¬ 
trum due to /nl), as a function of the multipole number 
1. The purpose of the figure is to illustrate how scale- 
dependent local non-Gaussianity can produce more than 
one signature in the CMB. Therefore, we have not in¬ 
cluded the effect of the Gaussian cosmic variance, which 
would add variance to both the power asymmetry and 
monopole modulation amplitudes; this becomes impor¬ 
tant when performing parameter estimation of /^l and 
n. While such a full parameter estimation analysis is 
beyond the scope of this work, we illustrate in the next 
section that, in a simplified context, adding the Aq con¬ 
straint can be useful in some situations. 

Notice that the scale-dependent non-Gaussian model 
of Eq. (45) also generates asymmetries in the spectral in¬ 
dex: 


dlnP 

dink 


ko 


{ns 


P)<f> + ■'T'/nl 


AP{ko,x) 

I -H AP(fco,x)’ 


(47) 


where {us — 1),^ is the spectral index of the Gaussian field 
and 


AP(fco, x) = 4 /°l I (48) 


is the super cosmic variance contribution. As before, AP 
may be expanded in multipole moments. A spatial mod¬ 
ulation of the spectral index was used in m as an al¬ 
ternative way to produce the power asymmetry. Here we 
see that such a modulation is a natural consequence of 
non-Gaussian scenarios that generate a scale-dependent 
power asymmetry. 
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FIG. 9. (color online). The effect of a simple scale-dependent 
local non-Gaussian model, /nl(^) = 50(I!/60)“° ®^, in the bis¬ 
pectrum (top), the power asymmetry amplitude (middle), and 
the modulation of power spectrum amplitude (bottom). A 
different choice of /Sl will only rescale the lines above. 


C. Statistical analysis 

The Planck 2013 analysis reported a power deficit 
at £ ^ 40 with a statistical significance of approxi¬ 
mately 2.5(7 — 3(7 [63], while the more recent Planck 
2015 analysis (SD] shows a slightly lower statistical signif¬ 
icance to the dehcit. Although that measurement alone 
does not require a new model, the scale-dependent non- 
Gaussian scenario generically generates a power modu¬ 
lation. Therefore, the data should be included in con¬ 
straining the model. If we include a measurement of the 
monopole modulation Aq in addition to the dipole power 
asymmetry amplitude A, the combined likelihood for /nl 
is: 


mNhlAA) =Tx(^’'^)Woid(^o,o-“™°). (49) 

In Figure [TO] we show examples in which, in addition 
to the power-asymmetry amplitude A, we also consider 
the monopole-modulation amplitude Aq. For simplic¬ 
ity, we assume that the result of a scale-dependent non- 
Gaussianity produces some large (statistically significant 
at « 3.3(7, for example) power asymmetry at large scales 
£ < 100 and which becomes smaller in magnitude and 
therefore less significant at small scales. In Figure [T^ 
therefore, we always take A = 0.055 and we have also 
added the estimate of /nl ~ —100 ± 100 at these scales. 


This simplification allows us to use the Gaussian cosmic 
variance from our constant /nl GMB realizations and 
simply augment the analysis in the previous section with 
the added information from a possible monopole mod¬ 
ulation. This is because we expect the scale-dependent 
non-Gaussianity to produce a monopole modulation Aq 
at large scales, the variance of which depends on the 
strength of non-Gaussianity at the scales used to obtain 
the monopole modulation amplitude Aq. We will define 
the measured Aq as in Eq.(36) i.e. a weighted average 
over multipoles 2 < £ < 100. Depending on the val¬ 
ues of ^/nlj ^ scale-dependent non-Gaussianity 

may generate significant monopole modulation at higher 
multipoles. But, we will restrict our analysis to a single 
value Aq obtained from the range 2 < £ < 100. Thus, our 
simple analysis is only sensitive to the average amplitude 
of non-Gaussianity at these scales and cannot constrain 
^/nl for which we will need to consider multiple bins of 
Aq in the data. 

For Aq [in Eq.(36l] generated at large scales from 
a scale-dependent non-Gaussianity, is the angu¬ 

lar power spectrum at large multipoles where the non- 
Gaussianity and therefore the effect of superhorizon 
modes on the observed power spectrum is small. In the 
case of the Planck satellite measurement of Aq in [63] , the 
best-fit angular power spectrum, which is also dominated 
by the larger-.^ modes, is taken as the 

In Figure]^ we see that for Aq — 0.02, which is about 
1.5 ct in the distribution of Aq values for our Gaussian 
GMB realizations, the posterior distribution shows more 
support for the Gaussian model. However, for a larger 
value Aq = 0.04 (about 2.9a), the addition of the Aq 
data favors the non-Gaussian model for a superhorizon 
efolds of iVextra = 50. For a larger number of super¬ 
horizon efolds Nextia = 100, the support for the non- 
Gaussian model again decreases because of the increased 
variance in the prediction from the non-Gaussian model 
itself. (However, some care should be taken in extrapo¬ 
lating this result since our analytic expressions are not 
valid when the non-Gaussianity in the large volume is 
strong.) The SDDR values for these cases are shown in 
Tabled 


-^0,obs 

-^extra 

SDDR (Boi) 

In Bqi 

- 

- 

0.2012 

-1.603 

0.02 

50 

1.107 

0.102 

0.04 

50 

0.0664 

-2.712 

0.04 

100 

0.1475 

-1.914 


TABLE II. SDDR for A — 0.055 and different observed val¬ 
ues of monopole power modulations Aq at large scales. See 
the text for a discussion on why the evidence increases and 
decreases for different values of Aq and Wxtra. 

We have only shown examples of how using Aq 
(i.e. power modulation from scale-dependent non- 
Gaussianity) can change the posterior probability (and 
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FIG. 10. (color online). The posterior distribution of |/nl| 
for Aobs = 0.055 and different assumed values of To.obs and 
A^extra. Note that the use of Ao changes the shape of the pos¬ 
terior distribution, and generates a peak at around /nl ~ 10 
for Ao = 0.04, A^extra = 50 (blue dotted line) in addition 
to the peak at |/nl| ~ 200 from the use of A = 0.055 
and the bispectrum constraints (black dashed line). For a 
larger A^extra = 100 (green solid line), the peak shifts towards 
|/nl| = 0 as the variance of the prediction of the non-Gaussian 
model increases; the evidence for the non-Gaussian model also 
decreases. 

hence the evidence) of scale-dependent non-Gaussianity 
in a simplified situation. A more involved analysis (i.e. 
including larger £ modes and measuring the scale depen¬ 
dence of /nl from bispectrum measurements with the 
latest data) is necessary to see how useful the inclusion 
of Aq, in addition to the power asymmetry amplitude, 
proves to be in the actual CMB data. In addition, we 
have chosen a simple model that may capture the features 
seen in the data but that may not be easily achievable 
from an inflationary point of view m- 


V. SUMMARY AND CONCLUSION 

We have performed a systematic study of the power 
asymmetry expected in the CMB if the primordial per¬ 
turbations are non-Gaussian and also exist, with no spe¬ 
cial features, on scales larger than we can observe. We 
have derived an expression to compute the expected de¬ 
viations from isotropy of the observed two-point function 
due to mode coupling with our specific realization of su¬ 
perhorizon modes, which generate anisotropy across our 
Hubble volume. Although we have focused our analy¬ 
sis on local non-Gaussianity, our method is quite general 
for describing deviations from statistical isotropy in a fi¬ 
nite sub-volume of an isotropic (but non-Gaussian) large 
volume. 

Exploiting the fact that local non-Gaussianity natu¬ 


rally produces a power asymmetry, we have shown how 
the observed asymmetry can be used for parameter esti¬ 
mation of the non-Gaussian amplitude /nl • We have also 
combined /nl constraints from bispectrum and power- 
asymmetry measurements to evaluate the Bayesian evi¬ 
dence for /nl ^ 0. In our simple examples, we find that 
the observed CMB power asymmetry only provides weak 
evidence for non-Gaussianity on large scales. 

Many previous works that propose mechanisms to gen¬ 
erate the CMB power asymmetry have required a de¬ 
parture from the simplest extrapolation of our observed 
statistics to larger scales in addition to non-Gaussianity. 
The most popular new feature is a large-amplitude super¬ 
horizon fluctuation that has an origin beyond inflation. 
Generating the observed scale dependence of the power 
asymmetry requires an additional elaboration, either in 
the form of scale-dependent non-Gaussianity or isocurva¬ 
ture fluctuations. However, our results show that scale- 
dependent non-Gaussianity of the local type is sufficient 
to generate a dipole power asymmetry at large scales 
without invoking exotic superhorizon fluctuations. A 
value of /nl consistent with the rather weak large-scale- 
only constraint (but well above the scale-independent 
bound) is enough to make the observed power asymme¬ 
try no longer anomalous. It is also worth noting that 
scale-dependent bispectra are only required for the local 
ansatz: scale-independent non-local bispectra can gen¬ 
erate a scale-dependent power asymmetry. All that is 
required to qualitatively match the observed scale de¬ 
pendence of the asymmetry is that smaller scales couple 
more weakly than large scales do to near-Hubble scale 
modes. 

For scale-dependent non-Gaussianity, we have demon¬ 
strated that a large-scale suppression or enhancement of 
the isotropic power spectrum is expected. Gonsequently, 
observations of a scale-dependent shift in the amplitude 
of the power spectrum can provide additional evidence 
in favor of non-Gaussianity. Although we have not found 
that the observed power deficit at large scales is suffi¬ 
cient to say the data strongly prefer the non-Gaussian 
model, we have also not performed a thorough analysis 
nor forecast the improvement possible with future data. 
Our parameter estimation and model comparison works 
were performed in a simplified setting to illustrate that 
a more careful look at the CMB data is necessary and 
useful. 

Our results strongly suggest that to evaluate the obser¬ 
vational evidence for models of the cosmological pertur¬ 
bations, departures from Gaussianity and from statisti¬ 
cal isotropy in the data can and should be considered to¬ 
gether. Furthermore, in weighing whether the large-scale 
CMB anomalies are evidence of something beyond slow- 
roll inflation lurking just outside the horizon, we should 
be sure to account for the possibility of non-Gaussian cos¬ 
mic variance contributions to anisotropy, which do not 
require the introduction of additional physics. 
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or equivalently, ((/)(k)(/)*(k')) = (27r)^5£)(k — \i')P^{k), 
and the dimensionless power spectrum V^j, is defined as 
V^{k) = 2TT'^P^{k)lk\ 

Appendix B: Derivation of modulated power 
spectrum for local non-Gaussianity 


For a homogeneous and isotropic cosmology, the two- 
point correlation function depends only on the magni¬ 
tude of the separation between the two points and is the 
Fourier transform of the power spectrum: 


(.(x-e.(x.I)> 


(fik 




(Bl) 


Appendix A: Gonventions and Definitions 

Throughout the paper, we follow the following Fourier 
convention: 


<('(x) = ^ 



(j){k) = ^ 


(Al) 


As usual, the temperature fluctuations of the CMB are 
decomposed into spherical harmonics: 

^ 

^ (fl) — ^ ^ (-^^) 

with the spherical harmonics normalization given by 

( 77 .) = Sl£/ 6 jn'm' • 

Therefore, the multipole coefficients of the CMB 
temperature fluctuations are given by 

f AT 

agm = J dnf,Yl^{h) — {h). (A4) 

For Sachs-Wolfe temperature fluctuations = 

— $/3) and when the Bardeen potential $ is a Gaussian 
field (d)(x) = (/>(x)), the Sachs-Wolfe aem is given by 

Att r d^k 

aem = —J j^(j){k)je{kx)Yl^{k). (A5) 


A hemispherical power asymmetry (a dipole modulation) 
cannot be described by allowing the power spectrum to 
be anisotropic in the usual way, P{k) —)■ P{k). (Since 
the fluctuations are real, the Fourier modes of k and —k 
are related, which forbids a dipole modulation. See, e.g., 
[82)1. Instead, if the amplitude of the coincident two- 
point function {(j)(x)(j){'x.)) varies spatially, we can likewise 
introduce position dependence into the power spectrum, 
P$(fc) —>■ P$(fc,x). Both even and odd multipole modu¬ 
lations can now be incorporated. 

This inhomogeneous power spectrum is well defined for 
Fourier modes with wavelengths that are much smaller 
than the length scale on which the power spectrum 
changes because these modes have unambiguous and con¬ 
stant wavelengths even if their amplitude varies spa¬ 
tially. In contrast, a plane wave with an amplitude 
that varies significantly within one wavelength doesn’t 
have a well-defined amplitude or wavelength. There¬ 
fore, Pj.(fc,x) is only defined for wavenumbers such that 
P$(/c,x) ~ P$(/c, x-f (27r/A:^)k). On a similar note, defin¬ 
ing 

(B2) 

only makes sense if P<s>{k,x) ~ P$(fc,x ± r); oth¬ 
erwise, it is futile to describe the two-point function 
(d) (x — r/2) $ (x -I- r/2)) in terms of a single power spec¬ 
trum. 


From this, we can obtain the Sachs-Wolfe angular 
power spectrum Ce defined as: duiSmm'Ci. 

cr = f f (A6) 

where we have defined the power spectrum P^{kl) as 

(</)(k)<^(k')) = {2TTfSD{k + k')P^{fc), (A7) 


In this appendix we show that a spatially varying 
power spectrum defined as in Eq.(B2) arises naturally 
in non-Gaussian scenarios where short-wavelength modes 
are coupled to long wavelength modes. For example, 
suppose the statistics in some large volume Vl ^ 00 
are described by the local ansatz in real space: $(x) = 
())(x) -k /nl (<(’(x)^ - ('(>(x)2)) where ij) is a Gaussian ran¬ 
dom field. Fourier modes of the non-Gaussian field are 
related to those of the Gaussian field by 
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$(k) = (/.(k) +/nl J ^^^[(/-(k-q)(/.(q)-(0(k-q)0(q))]. 


(B3) 


The two-point correlation of the non-Gaussian field is 

, - , , - , , . /■ d3k' 

$ X - - " 


(-I)K-I)) 


(27r)3 J (27r)3 


($(k)$(k'))e 


''i\pdk-l-k')-Xgi(k-k')'§ 


(B4) 


J (27r)3 

-h fm f f hk+k')-x i(k-k'). 

+ (2^)3 7 (27r)3" 


(27r)3 

+o{fL'Pl) • 




(27r)3 


When we consider statistics entirely within an infinite 
volume, the terms proportional to a single power of /nl 
vanish because is a Gaussian held. In that case, the 
power spectrum is corrected only by the last term (pro¬ 
portional to /nl)i which is small when non-Gaussianity 
is weak. 

Now, suppose we instead consider the statistics in a 
single sub-volume. Each sub-volume sits on top of a sin¬ 
gle realization of modes with wavelengths the size of the 
sub-volume or larger. So, for example, if mode p corre¬ 


sponds to a long wavelength mode while modes k and k' 
are well within the sub-volume, then 


((^(k )((>(k)(^(p)) |sub—volume — l/^(p)(0(k )0(k)) |sub—volume- 

(B5) 

That is, (/>(p) takes a particular value in the sub-volume, 
while ^(k) and ^(k') are still randomly distributed. 
Then, taking Eq. (B4) and considering a Eourier mode 
(/)(p) to be stochastic only if p = |p| > k^in- 




d^k 


sub—volume 


'|k|>fe„i„ (27r) 

(Pk 




(B6) 




bp|<fcmin 


+ o{fLPl), 
r Pk 

7|k|>fe„i„ (27r)3 




1 + 4/nl 


qp|<fc„ 


d^p 

(27r)3 


</'(p) cos (^) 


p^P X 


(B7) 


As discussed above, P<i,(fc,x) can be dehned in terms of 
the two-point function ($ (x — r/2) $ (x -|- r/2)) only if 
P$(fc, x) ~ P<!>{k, X ± r). In Eq.( |B7| ), we see that the spa¬ 
tial variation of the power spectrum arises from the 
factor, which implies that P$(fc,x) will be nearly con¬ 
stant on scales that are shorter than I/|p|, i.e. scales that 
are well within the sub-volume. Therefore, we should 
restrict the two-point function to separations such that 
p • r ^ 1, in which case the cosine term in Eq.(B7) is 
approximately unity. Then Eq. (B2) implies that 


P<i>,s{k,x) ~ P^{k) 


1 + Vni./ 


where the subscript £ specifies that the integral is only 


over long wavelength modes (fc < fcmin)- 

Expanding the factor we can see the effects of 

non-zero long wavelength modes as a multipole expan¬ 
sion: 


P<^,s(k) = P4k)[l + h^goo (B8) 

+/nl ^ 5'iM4iM(d-) + • • ■ ] • 

M={-1,0.1} 


The monopole shift is not observable, so we absorb it 
into the coefficient to define the observed isotropic power 
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spectrum: 

Pt^ik) = [1 


(B9) 


/nl \ 
1 + /nlSoo / 


giAfYiMin) 

M={-1,0.1} 


Finally, the shift to the power spectrum also shifts the 
observed value of /nl as defined from the local template 
for the bispectrum: 


($s(fci)$s(fc2)$s(A:3)) = {^TTfSoiki + k2 + k3)B°'^^{ki,k2, ks); (BIO) 

ks) = 2/nl[1 + hLgoo]P^{ki)P 4 ,ik 2 ) + sym, 

= 2f^l^Pt%k,)Pt^{k2)+sym. 


Since PI*"® = ^^[1 + /nl5oo] (considering only the 
isotropic piece), the second and third lines imply that 
/Sif = /nl/[1 + /nl5oo]- Then, the expression for the 
power spectrum, including the dipole asymmetry, is 


Related useful works that derive estimators when 
the primordial temperature field is modulated by an 
anisotropic field include [83l |84] . 


Ptik) = P; 


lObs 


(fc) 


1 +/nl* y]] giMYiM{n) + . ■ ■ 


M={-1,0.1} 

(Bll) 

The superscript “obs” on /nl should also be taken to 
indicate that the value does not contain the contribution 
corresponding to the Maldacena consistency relation for 
single clock inflation [/nl oc (n^ — 1)], which is unobserv¬ 
able [50] , 


Appendix C: Bipolar Spherical Harmonics 

Many works quantifying the likelihood of the power 
asymmetry make use of bipolar spherical harmonics 
[5E115S] . As an aid to that analysis, we repeat the calcula¬ 
tion of Appendix but for the atm ■ In the Sachs-Wolfe 
approximation, the statistics of the observed CMB two- 
point function are simply related to the two-point corre¬ 
lation of the potential at the time of decoupling: 


= IJ dn^YlmAni) J (^ 2 ) 




(27r) 


—ik2-ri2a: 


($(ki)$*(k2)), (Cl) 


where x is the comoving distance to the last scattering surface. When the potential is non-Gaussian according to the 
local ansatz, Eq.(B3), the statistics observed in a sub-volume will depend on the realization of the long wavelength 
modes: 


)sub—volume g ^ (^l)^2m2 (^2) 


qk|>fe„ 


d^ki 

(^ 




(C2) 


+2/nl 

+2/nl 


d^ki 


'|k2|>fcmin (271')^ 


,Piki) 

P{k2) 




/Pl<fer 

/ 

y |p|<fe, 


(27r)3 

iin ' ' 

d^p 

„i„ (27r)3 


^*(p)g*ki-nia:g-i(ki-|-p)-n2a: 


pi(p+k2)'nia;g—jk2'ri2a: 


= +2/NL(47r)(-l)-^ E 


Ci, 


'|p|<fe„ 


d^p 

(27r)3 


(j)*{p)YLM{p) + Q2 


qp|<fc„ 


d3p 

(27r)3 


(t>{p)YLM{p) 


Here Ce is given by Eq.(l30|) as usual, and are the Clebsch-Gordon coefficients. 
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The standard notation for the bipolar spherical harmonic expansion is 


LM 


so that the local model gives 


= 2/NL(4i)(-.)^/Shihei^c„'‘'-^ 


47r(2L + 1) 


-'000 


(C3) 


(C4) 




' |p|<fcn: 


d^p 

(27r)3 


(l)*ip)YLM{p) + Ce^ 


Appendix D: Notes on numerical realizations 


In Section III A[ we used Gaussian and non-Gaussian 
Sachs-Wolfe GMB maps to test some of our analytical 
formulas. The generation of non-Gaussian GMB maps for 
a more realistic GMB sky has been described in [HZIIHH]- 
Following |88j . we can write the temperature multipole 
moments as an integration over multipole moments of the 
Bardeen potential as a function of comoving distance r: 

atm = J r'^^em{r)at{r), (DI) 

at{r) = ^ j dk gi{k)jt{kr), (D2) 

where gg{k) is the transfer function of temperature in mo¬ 
mentum space. The process of generating non-Gaussian 
GMB maps as detailed in [SS], therefore, requires generat¬ 
ing non-Gaussian ^tm at different comoving distances 
considering the covariance ($£^^1 i'i^i)^i 2 m 2 (^ 2 ))) and nu¬ 
merically integrating over the comoving distance r. For 
a Sachs-Wolfe universe, gt{k) = —jt{krcnih)/5, and [8^ : 

2 f°° 

oii{r) = / dk k‘^jt{krcmh)jt{kr) 

371- Jo 

5D{r - Tcmb) 


Q„2 


Using this, we obtain: 


agm — 2 (7’cmb ) ■ 


(D3) 


Therefore, we can generate the local non-Gaussian 
GMB Sachs-Wolfe temperature anisotropies simply using 
Eq. (291. The maps used in our study were generated us¬ 
ing the HEALPIX software [90] , with £max = 300 for the 
input Cgs and = 128 for map-making. With the 

10000 Sachs-Wolfe simulated GMB maps for Gaussian 


bpi<fct. 




(I>(j>)Ylm{p) 


I 

and non-Gaussian potentials, we can perform a number 
of correlation tests in the Gaussian and non-Gaussian 
maps to understand what generates the power asymme¬ 
try. 

In Figure[^ we show three plots to illustrate some use¬ 
ful correlations among quantities in the simulated maps. 
In the top panel of Figure 11 we plot the directional 


power asymmetry amplitudes Ai for the Gaussian maps 
/nl = 0 and non-Gaussian maps with /nl = 500. The 
significant correlation between /nl = 0 and /nl = 500 
(with a correlation coefficient, p = 0.5) is expected be¬ 
cause both maps contain approximately the same amount 
of power asymmetry that comes from Gaussian cosmic 
variance. The /nl = 500 GMB skies contain additional 
power asymmetry on top of the /nl = 0 asymmetry. We 
have checked that the correlation gets weaker for larger 
values of /nl- 

In the middle panel, we plot Ci (which we have used 
to model the background dipole anisotropy in density 
fluctuations) against the dipole asymmetry amplitude A 
measured in both Gaussian (/nl = 0) and local non- 
Gaussian (/nl = 500) models. The Ci values remain 
almost the same with a correlation coefficient of 0.995; 
this indicates that the contribution from the 
term to the power spectrum is small for /nl = 500. We 
find that the correlation between Ci and A for the Gaus¬ 
sian maps is very weak (correlation coefficient = 0.035) 
compared to the correlation in the /nl = 500 model (cor¬ 
relation coefficient = 0.685). This shows that the com¬ 
bination of the background dipole anisotropy and non- 
Gaussianity is responsible for the power asymmetry in 
the non-Gaussian maps. To further test this notion, we 
also generated a set of /nl = 500 non-Gaussian maps in 
which Cl = 0. As we can see in the bottom panel of Fig¬ 
ure [m we do not find a significant increase in the power 
asymmetry distribution A for /nl = 500 with Ci = 0, 
which provides further evidence that the non-Gaussian 
power asymmetry is generated by a mode coupling be¬ 
tween a background dipole anisotropy and the small-scale 
modes. 
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FIG. 11. (color online). Top: The relation between the Ce 
asymmetry in three orthonormal directions for 1000 of our 
simulated Gaussian maps (x-axis) and the corresponding lo¬ 
cal non-Gaussian maps with /nl = 500 (y-axis). We find a 
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correlation between the power asymmetry amplitude and Ci 
for Gaussian maps (blue circles) and non-Gaussian maps with 
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try amplitude with Ci in the /nl = 500 realizations is strong. 
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